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We study the two-dimensional structure of thin transonic accretion discs in the 
vicinity of a non-spinning black hole within the framework of hydrodynamical ver- 
sion of the Grad-Shafranov equation. Our analysis focuses on the region inside 
the marginally stable orbit (MSO), r < r^^. We show that all components of the 
dynamical force in the disc become significant near the sonic surface and (espe- 
cially) in the supersonic region. Under certain conditions, the disc structure can 
be far from radial, and we review the affected disc properties, in particular the role 
of the critical condition at the sonic surface. Finally, we present a simple model 
aimed at explaining the quasi-periodical oscillations that have been observed in 
the infra-red and X-ray radiation of the Galactic Centre. 



1. Introduction 

The investigation of accretion flows near black holes (BHs) is undoubtedly 
of great astrophysical interest. Substantial energy release must take place 
near BHs, and general relativity effects, attributable to strong gravitational 
fields, must show up there. Depending on external conditions, both quasi- 
spherical and disc accretion flows can be realized. The structure of thin 
accretion discs has been the subject of many papers. Many results were 
included in textbooks [1, 2]. Lynden-Bell [3] was the first to point out that 
supermassive BHs surrounded by accretion discs could exist in galactic 
nuclei. Subsequently, a theory for such discs was developed that is now 
called the standard model, or the model of the a— disc [4, 5, 6]. 

1 
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Since then the standard disc thickness prescription has been widely used, 



where the disc thickness H is assumed to be determined by the balance of 
gravitational and accreting matter pressure forces with the dynamical force 
neglected. This relation was later used in the renowned approach where 
all quantities were averaged over the disc thickness [7] , with a lot of such 
one-dimensional models following [8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. As for 
the two-dimensional structure of accretion discs, it was investigated mostly 
only numerically and only for thick discs [9, 18, 19, 20]. 

Even though standard disc thickness prescription (1) and the averaging 
procedure are likely to be valid in the region of stable orbits r > [1], 
they in our view require a more serious analysis. It is the assumption 
that the transverse velocity V0 may always be neglected in thin accretion 
discs all the way up to horizon [21], i.e. that the disc thickness is always 
determined by (1), that is the most debatable [22]. This assumption is 
widely used, explicitly or implicitly, virtually in all papers devoted to thin 
accretion discs [23, 24]. 

There is a brief discussion on thin discs in this volume that covers theory 
shortcomings and a brief description of our approach (c/. Sec. named Thin 
disk in [22]). In this paper we briefly describe our study of subsonic and 
transonic regions of thin discs followed in Sec. 5 by the elaborate discussion 
of the supersonic flow. Finally, in Sec. 6 we develop a toy model for ex- 
plaining the observed quasi-periodical oscillations detected in the infra-red 
and X-ray observations of the GC. 

2. Basic equations 

We consider thin disc accretion on to a BH in the region where there are 
no stable circular orbits. The contribution of viscosity should no longer be 
significant here [22]. Hence we may assume that an ideal hydrodynamics 
approach is suitable well enough for describing the flow structure in this 
inner area of the accretion disc. Below, unless specifically stated, we con- 
sider the case of non-spinning BH, i.e. use the Schwarzschild metric, and 
use a system of units with c = G = 1. We measure radial distances in the 
units of M, the BH mass. 

In Boyer-Lindquist coordinates the Schwarzschild metric is [25] 



ds^ = -a^dt^ -I- Qikdx'-dx^, 



(2) 
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where 

= 1 - 2/r, Qrr = a"^, gee = r^, g^^ = zu^ = r"^ sin^ 9. (3) 

Wc reduce our discussion to the case of axisymmetric stationary flows. 
For an ideal flow there are three integrals of motion conserved along the 
streamlines, namely entropy, S, energy E = ^la'y, and ^;-component of 
angular momentum L = fizuu^, where fi = (p„j + P)/n (p„j is internal 
energy density, P = nT is pressure) is relativistic enthalpy. The relativistic 
Bernoulli equation u? = 7^ — m| — 1, where is the physical poloidal 
4- velocity component [22, 26, 27] , now becomes 

4 = —■ 4 

Below we use another angular variable Q = — 9 and for the sake 
of simplification we adopt the polytropic equation of state P = k{S)nF so 
that temperature and sound velocity can be written as [1] 

T = k{S)n^-^; cl = -k{S)rf-^. (5) 
3. Subsonic flow 

Following Sec. 1, we assume that the a-disc theory holds outside the MSO. 
We adopt the flow velocity components, which this theory yields on the 
MSO r = r^j^g,^ as the first three boundary conditions for our problem. For 
the sake of simplicity we consider the radial velocity, which is responsible 
for the inflow, to be constant at the surface r = r^-^^ and equal to uq and the 
toroidal velocity to be exactly equal to that of a free particle revolving at 
r = Tj^g.^ We also assume the speed of sound to be constant at the MSO, 
Cg = Co = const. Having introduced [27] — the Lagrange coordinate 
of streamlines at the MSO for Cg -C 1, i.e. non-relativistic temperature, 
we obtain from (4) and (5), 

4 = 4 + ^' + f^^{cl-c'.) + l{^L-&') + --- (6) 

The quantity 

^ ' a2 - a2 9^3 ' 



nearly parallel inflow with a small radial velocity Vr ~ aggCg/fK <C <S 1; = 
3rg, where rg = 2M is the gravitational radius of the BH of mass M. 
^FoT a free particle revolving at r = around a non-spinning BH we have [25] 
"v{^ms) = 1/^3, ao = a(r-ms) = \/2Ato = 7(»-ms) = \/4/3. 
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where Cq = Eq/^q and la ~ L/{jia cosGj^^g), is the poloidal four- velocity of 
a free particle having zero poloidal velocity at the MSO. 

In the extreme subsonic case, <^ Cg, the Grad-Shafranov hydro- 
dynamic equation is significantly simplified [22, 27]. The numerical re- 
sults are shown in Fig. 1. In the subsonic region, r* < r < r^g = 3rg, 
the disc thickness rapidly diminishes, and at the sonic surface we have 
-ff(r*) = uo/coH{r^g), so that you cannot neglect the dynamical force 
there. 

We stress that taking the dynamical force into account is indeed ex- 
tremely important. This is because, unlike zero-order standard disc thick- 
ness prescription (1), the Grad-Shafranov equation has second order deriva- 
tives, i.e. contains two additional degrees of freedom. This means that the 
critical condition only fixes one of these degrees of freedom {e.g. imposes 
some limitations on the form of the flow) rather than determines the angular 
momentum of the accreting matter [22, 27]. 

4. Transonic flow 

In order to verify our conclusions we consider the flow structure in the 
vicinity of the sonic surface in more detail. Since the smooth transonic flow 
is analytical at a singular point r = r*, 6 = [28], it is possible to express 
the quantities via a series of powers of h= {r — r*)/rt, and 9. Substituting 
these expansions into equations of motion we get a set of equations on the 
coefficients which allows us to reconstruct the flow structure in the vicinity 
of the sonic point (see Fig. 2), in particular 

where = [2/(r + 1)] >' ^ ' cq/uq gives the compression of stream- 
lines, ao = H{r^g)/ H{ri,), and rji ~ Uq^. Equation (8) yields shape of the 
sonic surface, u- = c^; it has the standard parabolic form h = (r-|-l)r?i6^/3. 

Since the transonic flow in the form of a nozzle (see Fig. 2) has longitu- 
dinal and transversal scales of one order of magnitude [28], near the sonic 
surface we have 5r\\ « 6r±, i.e. Sr^^ w H{r^). Hence for thin discs {i.e. for 
Co ^ 1) this longitudinal scale is always much smaller than the distance 
from the BH, Sr\\/r^, « H{r^,)/r^, <C 1. Only by taking the transversal 
velocity into account do we retain the small longitudinal scale 5r\\ <C rg. 
This scale is left out during the standard one-dimensional approach. 



2 2 



„2 _ 2 



1 - 2r,,h + -(r - 1) ^e^ + -(r + l)riiQ^ 

D Cq O 

1 + (r - 1) + i(r - 1) ^e^ - ^(r - i)(r + 



(8) 
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5. Supersonic flow 

Since the pressure gradient becomes insignificant in the supersonic region, 
the matter moves here along the trajectories of free particles. Neglecting 
the VgP term in the ^-component of relativistic Euler equation [29], we 
have [21] 

«».2(^ + ^^ + M^t.ne = 0. (9) 

Here, using the conservation law of angular momentum, can be easily 
expressed in terms of radius: = 2\/3/r. We also introduce dimensionless 
functions /(r) and g{r): 6/(r) = rUg, and g{r) = —auf > 0. Using (9) 
and the definitions above, we obtain an ordinary differential equation for 
/(r) which could be solved if we knew g{r): 

df f + 12 



dr r^g{r) 

From (4) we have m? w"^ as r ^ r^. On the other hand, ?i. 



(10) 



Co for r < r*. Therefore, the following approximation should be valid 
throughout the rg < r < r* region, g{r) « -\/(mwp"+"(acZ)^. 

Equation (10) governs the supersonic flow structure for the case of non- 
spinning BH. To get a better match with observations (c/. Sec. 6), we also 
consider a more general case of spinning BH, i.e. a Kerr BH with non-zero 
specific angular momentum a. After some calculation, equation (10) can 
be generalized to the Kerr metric with a strikingly simple form, 

d/ ^ P+a^{l~el)+ll ^ 
dr r'^g (r) ' 

where g (r) is a straightforward generalization of g{r) to the Kerr case; we 
omit it here due to space limitations. For the Schwarzschild BH (a = 0, 
eo = vWQ, Iq = 2\/3 [28]) equation (11) reduces back to (10). 
Integrating (11), we obtain 

tt" 



f (r) = K tan 



egiO 2 



(12) 



where k = y/a"^ {1 — Cq) + /g and tt/2 has been to a good accuracy substi- 
tuted for the integration constant arctan [/(r*)/\/3] .'^ 

The results of numerical calculations are presented in Fig. 1. In the 



■^For r just below ? », the function / should be positive to reflect the fact that the flow 
diverges. Then, / = corresponds to the point where the divergency finishes, and the 
flow starts to converge. 
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Figure 1. The structure of a thin accretion disc (actual 
scale) for cq = 10~^, uo = 10~^ after passing the MSO 
r = 3rg (a = 0, Schwarzschild case). As sufficient dissipa- 
tion can take place in the vicinity of the first node r = ri , 
we do not prolong the flow lines to the region r < ri . 




Figure 2. Schematics of 
thin disc streamlines pro- 
file around the sonic point. 
The flow has the form of 
the standard nozzle. Here 
X = ~h, and y = Q. 



supersonic region the flow performs transversal oscillations about the equa- 
torial plane, their frequency independent of their amplitude. We see as 
well that the maximum thickness of the disc in the supersonic (and, hence, 
ballistic) region, which is controlled by the transverse component of the 
gravitational force, actually coincides with the disc thickness within the 
stable orbits region, r > r^^, where standard estimate (1) is correct. 

Once diverged, the flow converges once again at a 'nodal' point closer 
to the BH. The radial positions the nodes are given by the implicit formula 
/(r„) = ±oo, i.e. 



(13) 



where n is the node number; the node with n — corresponds to the sonic 
surface. In this formula the sonic radius = tq can be to a good accuracy 
approximated by Tj^^ = '"ms('^) the expression for which can be found in 
most textbooks [1]. Figure 3 shows the positions of nodes for different 
values of cq (the positions do not depend on uq for uq <C cq) and the BH 
spin parameter a. The matter travel time between the nodes has weak 
dependence not only on uq but also on cq as well. This provides a means 
for testing the theory via observations, and we do this in the following 
section. 
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Figure 3. Radial positions (in the units of M) of the 
nodes for a range of initial sound velocities. Dotted, 
dashed, and solid curves correspond to the cases a = 
—0.5, a = 0, and a = 0.5 respectively. Each curve relates 
the radial position of a node to a value of the initial sound 
velocity. Intersection points of these curves with the line 
CO = const give the the nodes' radial positions for that 
particular value of cq . 



Figure 4. Relation of a vs. 
M. Dashed, dash-dotted, 
and solid lines come from 
matching T (700 ± 100 s), 
Tms (2200 ± 300 s), and To 
(1100 ± 100 s) respectively. 
The resulting error polygon 
is bolded. 



6. Applications to observations 

Suppose some perturbation in the disc (a "chunk") approaches the MSO. 
We expect to observe radiation coming from the chunk with the period of 
its orbital motion, 

T^s (a) = 27r (rfl^ + a) , (14) 

where a is the angular momentum per unit mass of the BH and r^^ is an 
estimate of the distance from the BH to the chunk [1]. After a number 
of rotations, the chunk reaches the MSO and passes through the nodal 
structure derived earlier (c/. Sec. 5) generating a flare. Each time the 
chunk passes through a node, it generates some additional radiation, and 
therefore the flare is likely to consist of several peaks. We believe that it is 
these peaks that were discovered in the infra-red and X-ray observations of 
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the GC [30, 31]. 

The time interval between the detection of two subsequent peaks equals 
the time it takes for the chunk to pass between two adjacent nodes (n-th 



where n is the index of the observed time interval (counting from one). 

The first term in r.h.s. of (15) can be easily obtained from the analysis 
of particle's geodesies in the equatorial plane [1] 



where the coefficients of the inverse metric are <?** = — S^/(/9^A) and i?**^ = 
wgi**; the definitions of S, p, A, and w can be found elsewhere in this 
volume [22]. 

For definiteness and simplicity, we assume that the observer is located 
along the rotation axis of the BH. On its way to the observer, the radiation 
travels along the null geodesic that originates at a node in the equatorial 
plane (e.g. r = r-n, = 7r/2) and reaches the observer at infinity [r = oo, 
9 = 0). Using these as boundary conditions for null geodesies in the Kerr 
metric [32], we numerically find T„ '{a^c^). 

Figure 5 shows the dependence of observed time intervals on the value of 
the speed of sound in the disc. Although each individual time interval may 
depend on cq, the range [Tniin(a), rmax(a)] of observed time intervals (see 
the caption to Fig. 5) is independent of cq. With such weak dependence on 
the speed of sound in the disc, we have only two matching parameters: the 
specific spin a and the mass M of the BH. 

In the flare precursor section we associate the period Tms with the 
2200 ± 300 s one (group 5, c/. Table 2 in [30]) and the time interval Ti 
with the period of 1100 ± 100 s (group 4 in [30]). In consistency with the 
infra-red observations of the flare, the periods Ti, T2, etc. chirp with the 
peak number [31], i.e. resemble the QPO structure and thus form a cumula- 
tive peak of a larger width shifted to higher frequencies on the flares' power 
density spectra (700± 100 s, group 3, c/. Fig. 3a and 4a in [30]). We can es- 
timate the average frequency of this peak as 1/T = 1/2 (l/Tmin -I- l/Tmax)- 
The results of the periods' matching procedure are shown in Fig. 4. Despite 
large uncertainties in the observational data allowing significant freedom of 
a and M, high positive values of a {i.e. the disc orbiting in the same direc- 
tion as the BH spin) are clearly ruled out. 





(16) 
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Figure 5. The dependence of time intervals between the peaks in a flare on the speed 
of sound in the disc, cq, for a moderately spinning BH (a = —0.5). The uppermost 
curve corresponds to the time interval Ti between the the 0th and 1st peaks, the second 
curve from the top corresponds to the time interval T2 between the 1st and 2nd ones, 
etc. All intervals T„ behave very similarly: they first decrease with cq and then abruptly 
increase to infinity due to time dilation when the innermost node in the pair comes close 
to the BH horizon (which is indicated with upward arrows). Even though individual time 
intervals between subsequent peaks in the flare may depend on the temperature in the 
disc (which is proportional to Cq, see (5)), their minimum and maximum values remain 
the same for the range of sound velocities where there are several intervals observed. 
In the particular case of a = —0.5, illustrated in the figure, we have Tmin ~ 32 and 
Tinax ~ 60 with all other time intervals lying in between. 
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